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Abstract 



a,: 
a: 

We discuss the polarised fragmentation functions of quarks and 
gluons in Perturbative Quantum Chromodynamics. The Altarelli- 
Parisi evolution equations governing these fragmentation functions are 
^ ■ presented. We find that the first moment of the polarised gluon frag- 

. mentation function to order a s is as important as that of the quark 

and anti-quark fragmentation functions. The A production in po- 
larised e + e~ annihilation where this can be realised is discussed. We 
propose the factorisation method to compute the gluonic contribution 
to A production to order a s . With appropriate operator definitions 
for the quark and the gluon fragmentation functions, we find that the 
hard part of the factorisation formula is free of any infrared singular- 
ities confirming the factorisation. 



Quantum Chromodynamics(QCD) is the most successful theory of strong 
interaction physics describing the dynamics of the quarks and gluons inside 
the hadrons |T|. Parton Model description of various strong interaction pro- 
cesses has been a successful description. Both inclusive processes such as 
Deep Inelastic Scattering(DIS) and hadro production in e + e~ annihilation 
can well be described in the QCD improved Parton model. In DIS, the 
probability distribution functions of quarks and gluons inside the incoming 
hadrons are the non-perturbative inputs of the model. Though they are non- 
perturbative, their evolution in terms of the scale is completely described by 
the well known Altarelli-Parisi(AP) evolution equations. In the processes of 
hadroproduction in e + e~ annihilation, the fragmentation functions of quarks 
and gluons into hadrons are the non-perturbative inputs but again their evo- 
lution is completely governed by the AP evolution equations. 

Needless to say that Perturbative QCD(PQCD) improved Parton model 
unravelled both momentum and spin structure of the proton in terms of its 
constituents such as quarks and gluons. The momentum structure of the 
proton was studied in the unpolarised lepton proton DIS experiments [[L]]. 
Few years ago a series of polarised lepton proton DIS experiments unrav- 
elled the contributions coming from quarks and gluons to the proton spin 
[§]. In both the unpolarised and the polarised DIS experiments, a remark- 
able fact emerged which was that the gluons share significant amount of the 
momentum and the spin of the proton contradicting the naive expectation 
that most of the momentum and the spin are carried by the valence quarks. 
These experimental facts are well understood by using the AP evolution 
equations of both unpolarised and polarised quark and gluon distribution 
functions. These inclusive experiments are excellent tests of QCD improved 
Parton model and the AP evolution equations. 

Recently, there has been a lot of works on the fragmentation functions of 
quarks and gluons measurable in the process of hadroproduction in e + e~ 
annihilation and also on the test of QCD improved Parton model subli- 
mented with the AP evolution equations for the fragmentation functions. 
All these analysises are on unpolarised e + e~ scattering experiments which 
require only the AP equations for unpolarised fragmentation functions. The 
QCD improved analysis of the polarised e + e~ scattering experiment is yet to 
be done. 

It is in this spirit, this letter presents a systematic analysis of the polarised 
fragmentation functions of quarks and gluons within PQCD. This analysis is 
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important as we will see that polarised gluonic contribution to the produc- 
tion of polarised hadrons is as important as polarised quark contribution to 
polarised hadrons. Our analysis is quite general in the sense that it can be 
applied to other inclusive hadroproduction processes such as polarised pro- 
ton proton and polarised proton anti-proton scattering experiments as these 
fragmentation functions and their evolutions are universal. 

Let D^q^(x, Q 2 ) be the probability distribution for a parton a with polar- 
isation h produced at a scale Q 2 to fragment into a hadron(iJ) of polarisation 
s carrying a momentum fraction x of the quark momentum. Now we define 
polarised parton fragmentation functions as 



If the strong interaction world consists only of quarks and a single hadron, 
then the fragmentation function is 5(1 — x) implying that all the spin will 
be transferred to the produced hadron. This is not usually the case because 
of two reasons. Firstly the spin of the parent parton can be distributed as it 
can undergo collisions which can flip the spin(mass corrections). Secondly, 
the single quark can fragment into many hadrons as it can excite the vacuum 
producing qq pairs. So the quark spin will be distributed among the produced 
hadrons. So it is not necessary that the quark spin is completely transferred 
to a single hadron. The situation is similar to DIS where the proton spin is 
shared by various partons. 

The evolution of the quark and gluon fragmentation functions are gov- 
erned by the AP evolution equations and are given by 



^ADjfof) = ^ £^ [AD»(y,t)AP qq (x/y) + ADf(y,t)AP gq (x/y) 

(2) 



where t = log(Q 2 /A 2 ) and a s (t) is the strong coupling constant. Here, 
a s (t)AP a b(y)dt is the probability density of finding a parton of type a at 
the scale t + dt with momentum fraction y inside the parton of type b at a 
scale t. The splitting functions are given by 



AD»(x, Q 2 ) = D™(x, Q 2 ) - D%»(x, Q 2 ) 



(1) 
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Here, C 2 (#) = (A^ 2 - l)/2iV, C 2 (G) = N and T(i?) = f /2 with N = 3 
for SU(3) and / being the number of flavours |L|. In the above equations, 
usual + prescription has been used. Notice that the splitting function matrix 
appearing in this AP evolution equation is just transpose of that appearing 
in the AP evolution equation for the parton probability distributions. This 
can be understood very easily: The emission of a quark(gluon) from a quark 
(gluon) only affects the probability of quark (gluon) fragmenting into hadron. 
On the other hand, the emission of a quark from a gluon changes the prob- 
ability of the gluon fragmenting into hadron. Similarly, the emission of a 
gluon from a quark affects the probability of quark fragmenting into hadron. 
That is why this matrix is transpose. These equations can easily be solved 
in the Mellin space. Define 

AD?(n,t) = / ar^ADf (x,t)dx 
Jo 

AP ab (n) = j\ n - l /\P ah (x)dx (4) 

It is interesting to consider the lowest moment n = 1. From eqn.(||), we find 
that the first moment of the polarised gluon fragmentation function to order 
a s (t) satisfies a simple first order differential equation, that is 

j^Df (t) = a.(t)A,A£f (f) (5) 

where @ = (11C 2 (G) —4T(R))/12tt. The solution to the above equation can 
be found very easily using renormalisation group(RG) equation for the QCD 
coupling constant, 

-a s (t) = -(3 a 2 s (t) (6) 
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From eqns.(P, ||), we obtain an interesting behaviour of first moment of gluon 
fragmentation function: the product of the first moment of polarised gluon 
fragmentation function times the strong coupling constant is scale indepen- 
dent to order a 2 s {t), 

j t (a s (t)ADf(t)) = 0(a s (tf) (7) 
In other words, to order a^(t), AD? increases as the scale t increases, i.e 

AD? (t)=K log (^j (8) 

where K is some constant. Recall that the counter part of such a relation in 
the polarised parton distribution functions exists and has opened up a better 
understanding of the spin structure of the proton ||. That is, due to the 
similar relation 

j t (a s (t)Ag(t)) = 0(a s (tr) (9) 

where Ag(t) is the first moment of polarised gluon distribution function, the 
polarised gluonic contribution to proton spin is significant at very high ener- 
gies. This was experimentally observed and the result was that the gluon car- 
ries significant amount of the spin of the proton. Hence the relation(eqn.([7|)) 
suggests that at very high energies the polarised gluon fragmentation into 
polarised hadrons is as significant as polarised quark fragmentation into po- 
larised hadrons. So care should be taken when one interprets the data on 
polarised hadron production in polarised e + e~ annihilation processes. 

Now let us consider the first moment of the polarised quark fragmentation 
function into polarised hadron. From eqn.(fj), it turns out 

jAD* (t) = \a s {t)AD? (t) (10) 

From eqns. (|7[) and (|TI]), we find that AD?_{t) grows as t. So this is a re- 
markable behaviour. This is very different from first moment of the polarised 
quarks in polarised hadron where it is scale independent. 

It is instructive to find out the consequences of these results from theoret- 
ical point of view. Let us consider the polarised inclusive hadroproduction of 
e + e~ scattering. In the quark sector, the first moment of the cross section is 
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proportional to AD^(t). The eqn. ([TOD suggests that the leading behaviour 
of ADg(t) is logarithmic. In the gluonic sector, the first moment of the 
gluonic contribution is proportional to AD^(t). The coefficient function is 
proportional to log{Q 2 ) times strong coupling constant a s {t). So, from the 
eqn. (|7D, the first moment of the gluonic contribution is of the same order as 
that of leading order quark contribution. Hence the naive expectation that 
polarised gluonic contribution to polarised hadroproduction is next to lead- 
ing order effect is not correct. The above observation is very general as it is 
based on the AP evolution equation for polarised parton fragmentation func- 
tions. So the above analysis is applicable to other polarised hadroproduction 
processes. 

In the following we consider the polarised hadroproduction in polarised 
e + e~ annihilation [|J. Since the gluonic contribution is as important as quark 
contribution, we compute the polarised gluonic contribution to polarised A 
production. We work in the energy range where only process where 'photon' 
fragmenting into hadron is dominant. The reason why we consider inclusive 
polarised A production is that the cross section is large compared to that 
of other particles. Also it is easy to measure the polarisation of A ||. Our 
analysis can be extended to other hadro productions. The complete analysis 
including Z exchange channel is reserved for future publication ||. 

The inclusive polarised A production rate factorises as: 

da(s, Sl ) = -^L^( qi , q2 , Sl )(^)A7rW^p,s)-%- (11) 

where L^ u (qi, q 2 , Si) is the leptonic part arising from e + e~ annihilation into 
a photon of virtuality Q 2 and W^ u (p,q, s) is photon fragmentation tensor. 
The arguments of these tensors are momenta described in the figures 1 and 
2 and s, Si are the spins of A and electron respectively. The photon frag- 
mentation tensor contains all the information about the polarised quark and 
gluon fragmentation into polarised A. 

The operator definition of W^(q,p, s) is found to be 

W*(q,p,s) = ±J d 4 ^ < 0\J,(0)\A(p,s)X >< A(p, S )X\J v (£)\0 > 

(12) 

where J M (0 is the electromagnetic current, X is unobserved hadrons (where 
the summation over all X is assumed), q is the virtual photon momentum 
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given by q = qi + q2, p and s are the momentum and spin of the polarised A 
detected. 

The polarised A production cross section gets contribution only from 
antisymmetric parts of leptonic and photonic tensor. The antisymmetric 
part of the leptonic tensor is found to be 

L^(qi, <?2, si) = -2ie 2 s 1 e^ wal3 q^q^ (13) 

where again qi, q 2 are the momenta of the incoming leptons and s\ is the spin 
of the polarised lepton (electron). The photonic fragmentation tensor is not 
calculable in Perturbative QCD(PQCD) as we do not know how to compute 
the matrix element of em(electro magnetic) current between A states and 
the vacuum. But this can be parametrised using Lorentz covariance, gauge 
invariance, Hermiticity, and parity invariance. Hence the antisymmetric part 
of this photon fragmentation tensor takes the following form 0: 

W*(q,p, s) = — e^gV^z, Q 2 ) + —e^ Xu q x { s CT - ^-pA g£(x, Q 2 ) 
* p.q p.q \ p.q ) 

(14) 

where x = 2p.q/Q 2 , Q 2 = q 2 and s 2 = —1. Here the polarised fragmentation 
functions g^(x, Q 2 ) are real and Lorentz invariant, hence they are functions 
of x and Q 2 . 

We consider the following asymmetric cross section such that only the 
§i(x, Q 2 ) structure function is projected out: 

<Mtt - ti) 2 7T aA 2 

dxdCos(e) =<*qS*9A*,Q)Co8(9) (15) 

Where a = e 2 /4n, 9 is the angle between produced hadron and the incom- 
ing electron. Here, |f means that both incoming electron and the produced 
hadron are paralelly polarised and ]{ means that they are polarised anti- 
paralelly. From the above asymmetry it is clear that the polarised fragmen- 
tation function (Ji(x, Q 2 ) can be measured only through angular distribution 
of polarised A [|J. 

Recall that in the DIS, the hadronic tensor appearing in the cross section 
formula is factorised into a perturbatively calculable hard part and non- 
perturbative matrix elements of quark and gluon operators sandwiched be- 
tween hadronic states(called soft part) Pfl. This factorisation ensures that 
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the hard part contains no IR singularities. Following similar factorisation 
procedure, one can write the photonic fragmentation tensor in terms of com- 
pletely calculable hard part denoted by H^ u (y, Q 2 ) and non-perturbative op- 
erator matrix elements denoted by D^(z, Q 2 ) whose operator definitions will 
be given later. Hence, the factorisation formula for photon fragmentation 
function reads as 

W$'\x, Q 2 ) = £ I' ^Wg? (*/!/. Q 2 )^(y, Q 2 ) (16) 

a (h) Jx y 

where a runs over all the partons such as quarks, anti-quarks and gluons, 
h is the helicity of the parton. H®( h \z,Q 2 ) is the hard part of the parton 
differential cross section for the production of parton of type a with polar- 
isation h and energy fraction z = 2p.q/Q 2 . The soft part AD^$(y,Q 2 ) is 
the usual polarised parton fragmentation function whose operator definition 
and the physical interpretation in terms of the partons will be discussed in 
the following. The above factorised formula ensures that the H® u (z,Q 2 ) is 
free of any Infrared(IR) singularities. In other words, if H^ u (z, Q 2 ) is free of 
any IR singularities, the tensor W£,(x, Q 2 ) is said to be factorisable. In the 
following we shall show that such a factorisation exists and that it works to 
order a s in the gluonic sector. This is the by product of our analysis on the 
evaluation of first moment of polarised gluonic contribution to polarised A 
production. 

Let us first work out the operator definition for D£(x, Q 2 ) to leading 
order. As we know the photon fragmentation tensor W^ u (x, Q 2 ) is not com- 
putable in the perturbation theory due to the non-perturbative nature of 
hadrons(here A). The physical interpretation of the photon fragmentation 
tensor in terms of the parton fragmentation functions can be easily got if we 
work on the light cone || since we are interested in the limit Q 2 — > oo(see 
below). To make the computation simple and the physics more transparent, 
let us assume that the currents appearing in the photon fragmentation ten- 
sor are free field currents so that formal manipulations can be done with out 
explicitly evaluating them between the A states and the vacuum. Note that 
dominant contribution to this integral in W^ u (x, Q 2 ) (eqn.(fl"2|)) in the limit 
Q 2 — > oo comes from the light cone region £ 2 — > 0, otherwise the exponential 
factor will kill the integral in this limit. In this lightcone limit £ 2 — > 0, let 
us assume that only one of the field operators (say quark) of each current is 
responsible for the fragmentation whatso ever be the mechanism. Hence the 
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rest of the operator matrix elements can be computed between free parton 
(say anti-quark) states to the leading order. After a little algebra we find 
that the quark contribution to the structure function g A (x, Q 2 ) to leading 
order on the light cone turns out to be 

gf(x,Q 2 ) = 3-J2e 2 q AD° q A (x,Q 2 ) (17) 

X g 

where 

AD° q A (x,Q 2 ) = ±J de-e-*-* + /* < 0| 7+ 75^(0)|AX >< AX$(r)|0 > 

(18) 

where average over colour is implicit in the eqn.(^). Similarly one can 
find out the anti-quark contribution to g A (x,Q 2 ). Following the procedure 
adopted in reference |10[, we find that the AD q A (x,Q 2 ) measures the num- 
ber of up(T) polarised quarks fragmenting into up(f) polarised A minus the 
number of down (J,) polarised quarks fragmenting into up(f) polarised A. 

Generalising the above idea, one can express polarised quark, anti-quark 
and gluon fragmentation functions in terms of quark and gluon operators as 



follows 101, fill: 



AD A +q -(x, Q 2 ) = Z-] e-^ p+ ' x (< 0| 7+ 75^(0)|AX >< AX|^(D|0 > 
+ < 0|7 + 75^(r)|AX >< AA>(0)|0 >) (19) 



AD A (x,Q 2 ) = -£^ J e"^ + /* < 0\Gt a (0)\AX >< AX|G a + a (r)|0 > 

(20) 

where, sum over X is assumed and the usual field strength tensor. Also, 
averge over colour is implicit. We have also dropped the Wilson link operator 
as it becomes identity in the lightcone gauge A + = which we choose to 
work with. The A carries the momentum p and spin s. It is clear from the 
above operator definitions of the parton fragmentation functions that the 
matrix elements of these operators for hadronic states are not calculable. But, 
interestingly, they are calculable for the partonic states such as quarks, anti- 
quarks and gluons in the perturbation theory as one knows how these partonic 
operators act on the partonic states. This fact that they are calculable for 
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the partonic states is exploited to compute the hard part (H®(z, Q 2 )) of the 
factorisation formula. 

First we compute the leading order contribution coming from the oper- 
ators to the polarised fragmentation function g^(x, Q 2 ). Using the factori- 
sation formula and the normalisation that AD^(x, Q 2 ) — 8(1 — x) to lowest 
order, we find that the quark and the anti-quark contributions to g\(x,Q 2 ) 
come form Fig. 1. which we denote by g®(x,Q 2 ), 

g^(x, Q 2 ) = 3- £ e\ (a<(x, Q 2 ) + AD*(x, Q 2 )) (21) 
i 

From the above equation, it is clear that the polarised fragmentation 
function gi(x, Q 2 ) to leading order a s (Q 2 ) measures the polarised quark and 
anti-quark fragmentation functions weighted by appropriate charge square 
factors(e^). 

The a s (Q 2 ) corrections to g^(x,Q 2 ) come from two sources: 1. gluon 
bremstalung and virtual corrections to polarised quark and anti-quark frag- 
menting into polarised A, 2. the polarised gluon fragmenting into polarised 
A. Here, we are interested only in the evaluation of polarised gluonic con- 
tribution to the production. Using the factorisation formulafeqn. (|16|) , the 
gluonic contribution to gi(x, Q 2 ) can be formally written as 

Q 2 ) = -H 4 t ~n g (x/y, Q 2 )ADg(y, Q 2 ) (22) 
x q Jx y 

We again use the factorisation formula to compute H^ u (z,Q 2 ). H^ u (z,Q 2 ) 
gets contribution from order a s (Q 2 ) onwards. To compute this, we replace 
polarised A by polarised gluon in the eqn. (p~6|) . With this replacement 
we find that the left hand side of the eqn. flTBp is just the cross section 
for photon decaying into polarised gluon and a quark-antiquark pair. On 
the other hand, the right side involves the evaluation of quark operator for 
polarised gluonic states to order a s (Q 2 ) and gi q (z,Q 2 ) and AD^(z,Q 2 ) to 
lowest order (a s (Q 2 ) ). Hence, 

H g (z, Q 2 ) = W g (z, Q 2 ) - 3ADf Q 2 ) (23) 

Let us now compute W g (z,Q 2 ) appearing in the eqn.(p3|). W g (z,Q 2 ) can 
be computed from the polarised cross section for e + e~ annihilating into the 
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polarised gluon: 

dcr 9 (s, Si 



is 



1 



Q 



where, z = 2p g .q/Q 2 , p g is the momentum of the polarised gluon and 9 
is the angle between produced gluon and the incoming electron. H g (p g ,q) 
appearing in the above equation (eqn. ([24])) is found to be 

Q 



32(2 



7T 



dx 2 V^\M 9 \ 



2 
lii' 



(25) 



where x 2 = 2p q .q/Q 2 and the projector V^ v = ie^ u x a p g q a /2p g .q. The pro- 
jected matrix element square V g .\M 9 \ 2 is computed from the Fig. 2 (with 
gluon polarised) and is given by 



pg u \M 9 \ 2 



[IV 



3C 2 (R) 



A(27r) 3 e 2 aa s { m 2 Q 2 - st" 



7T V (s + t) 2 
Am 2 Q 2 - 2m 2 g s - 2Q 2 s + s 2 - t 
t 2 

Am 2 g Q 2 - 2m 2 g t - 2Q 2 t + t 2 -s 2 



{s + t)(s + t-m 2 g -Q 2 
si 



(26) 



where the Mandelstam variables s = (pq+p g ) 2 ,t = {p q +p g ) 2 and u = {p q +Pq) 2 
satisfy s + t + u = m 2 + Q 2 . Also note that s = Q 2 (l — x 2 ). Substituting 



the eqns. Q2HD and (^) in the eqn. (£23|), we find that the asymmetry in the 
polarised gluon emission turns out to be 



<Mtt - TP = 3C 2 (R)e 2 a 2 a 
dzdCos{9) 



Q 2 



(2 -z) log 



2(2 - z)(3 z 



Cos(9) 
(27) 

where z = 2p g .q/Q 2 , (3 Z = (1 — Am 2 /Q 2 z 2 ) 1 ^ 2 and 9 here is the angle between 
the outgoing polarised gluon and incoming lepton (electron). The above result 
shows that there are no soft singularities. The small nonzero gluon mass is 
used to regulate collinear divergence. Since we are looking at the polarised 
gluon production, there is no virtual correction to this order, hence there is 
no UV divergences. From the eqn. (27) , we find that 

,Q 2 z 2 , 



W g (z,Q 2 ) = 3C 2 (R)-^2 



2tt 



(2-z)log( 



mi 



2(2 



(28) 
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Now, let us compute the matrix element D^(z, Q 2 ) appearing in the eqn. 
(|23|). We can compute this using the eqn. (|19|) with A replaced by polarised 
gluon. The diagrams contributing to this matrix element are given in Fig.3. 
This is found to be 



AD 9 q+ g(z,Q 2 



AC 2 (R)a s z I -^%2ir8(k 2 )2n5(k + - p+ + p+/z) (29) 



Tr 



7 75 



(2ir) d ' 

iPg 



J 9 1 fg 
(k - pg) 2 + ie 



The above expression is UV singular and we regulate this divergence in Di- 
mensional Regularisation method. A small gluon mass m g is introduced to 
regulate IR divergence. From the detla functions appearing in the above 
equation we find that k + = p^(z — l)/z and k~ = k\z/{2p+(z — 1)). Using 
these constraints, the k + and k~ integrations can be safely done. The k\ 
integration can also be done and the result is 



[z,Q 2 )=C 2 (R)^2 



(2-z) log 



mg 2 (l 



(30) 



The UV divergence appearing in the eqn. (|30|) is renormalised in MS scheme, 
hence the appearance of substration scale ix 2 R . 

Substituting the eqn. (|30|) in eqn. (|2^) , we find that the gluonic coefficient 
function is 



H 9 {z,Q 2 ) 



3C 2 (R)—2 
V J 2n 



(2-z)log 




+ (2 - z) log (1 - z) + 3z - 4 



(31) 

where hr is the scale at which the gluon fragmentation operator matrix ele- 
ment is renormalised. From the above equation it is clear that the collinear 
divergence appearing in the cross section formula cancels against that ap- 
pearing in the matrix element leaving the hard part H g (z, Q 2 ) free of any IR 
divergences. This is the proof of factorisation in the gluonic sector to order 
a s (Q 2 ). Substituting the above equation in eqn. ([22]), we get the gluonic con- 
tribution to Qi (x, Q 2 ) and hence the gluonic contribution to the asymmetry 
given in the eqn. (|15|). 

In this paper we have systematically analysed the importance of gluons to 
the production of polarised A in an inclusive e + e~ annihilation. We have done 
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this using the AP equation for the gluon fragmentation function. We have 
computed the gluonic contribution to order a s to the polarised fragmentation 
function g± (x, Q 2 ) appearing in the asymmetry of the polarised cross section 
using the factorisation formula. The factorisation that the cross section can 
be factored into a hard and a soft part has been demonstrated in the gluonic 
sector. 
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Figure Captions: 

1. Graph contributing to e~(gi)e + (g 2 ) — > q(Pq)q(p q )- 

2. Graphs contributing to e~(gi)e + (g 2 ) — ► q(Pq)q(Pq)g(Pg)- 

3. Graph contributing to AD^(x, Q 2 ) 
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